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ABSTRACT 

Based on a large Fermi blazar sample, the blazar sequence (synchrotron peak fre¬ 
quency iZpeak versus synchrotron peak luminosity Lpeak) is revisited. It is found that 
there is significant anti-correlation between z/peak and Fpeak for blazars. However, after 
being Doppler corrected, the anti-correlation disappears. The jet cavity power (Cjet) 
is estimated from extended radio luminosity. So it is free of beaming effect. We find 
that there are significant anti-correlations between Pjet and beam-corrected zZpeak fo^' 
both blazars and radio galaxies, which supports the blazar sequence and unification 
of blazars and radio galaxies (an alternative relationship is the correlation between jet 
power and 7-ray photon index). 

Key words: BL Lacertae objects: general - galaxies: jets - quasars: general - gamma- 
rays: galaxies 


1 INTRODUCTION 

Blazars are the most extreme and powerful active galactic 
nuclei (AGN), pointing their jets in the direction of the ob¬ 
server (Urry & Padovani 1995; Ghisellini & Tavecchio 2015). 
Based on the equivalent width (EW) of the emission lines, 
the blazars are often divided into two subclasses of BL Lac¬ 
ertae objects (BL Lacs; rest frame EW< 5 A) and flat spec¬ 
trum radio quasars (FSRQs; Urry & Padovani 1995). Ghis¬ 
ellini et al. (2009, 2011) introduced a physical distinction be¬ 
tween the two classes of blazars, based on the luminosity of 
the broad line region measured in Eddington units. Giommi 
et al. (2012, 2013) suggested that blazars should be divided 
into high and low ionization sources. Their broadband spec¬ 
tral energy distributions (SEDs) are usually bimodal. The 
lower bump peaks at IR-optical-UV band and the higher 
bump at GeV-TeV gamma-ray band (Ghisellini et al. 1998; 
Abdo et al. 2010b). 

Fossati et al. (1998) presented the blazar sequence: the 
peak luminosity of the synchrotron component (Lpeak) and 
the Compton dominance (the ratio of the peak of the Comp¬ 
ton to the synchrotron peak luminosity) are anti-correlation 
with the synchrotron peak frequency Z/'peak. Ghisellini et al. 
(1998) modeled the broadband spectra of blazars, and sug¬ 
gested that along the blazar sequence (BL Lacs —^ FSRQs), 
a stronger radiative cooling causes a particle energy distri¬ 
bution with a break at lower energies. So the synchrotron 
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and Compton peaks shift to lower frequencies. Various lines 
of evidence support the blazar sequence (Georganopoulos 
et al. 2001; Cavaliere & D’Elia 2002; Maraschi & Tavec¬ 
chio 2003; Xie et al. 2007; Maraschi et al. 2008; Ghisellini 
& Tavecchio 2008; Ghisellini et al. 2009, 2010; Abdo et al. 
2010a; Chen & Bai 2011; Sambruna et al. 2010; Finke et al. 
2013; Chen 2014). But many contrary opinions have been 
presented (Giommi, Menna & Padovani 1999; Caccianiga & 
Marcha 2004; Padovani et al. 2003; Anton & Browne 2005; 
Nieppola et al. 2006, 2008; Padovani 2007; Giommi et al. 
2012). The authors mainly focus on sample selection effect 
and beaming effect. Ghisellini & Tavecchio (2008) revisited 
the blazar sequence and proposed that the SED of blazar 
emission are linked to the mass of black hole and the ac¬ 
cretion rate. Finke (2013) studied a sample of blazars from 
2LAC (Ackermann et al. 2011) and found that a correlation 
exists between Compton dominance and the peak frequency 
of the synchrotron component for all blazars, including ones 
with unknown redshift. And via constructing a model, Finke 
(2013) reproduced the trends of the blazar sequence. Meyer 
et al. (2011) revisited the blazar sequence and proposed the 
blazar envelope: FR I radio galaxies (FR Is) and most BL 
Lacs belong to weak jet population while low synchrotron 
peaking blazars and FR II radio galaxies (FR IIs) are strong 
jet population. 

Many evidences support that FR IIs are the parent pop¬ 
ulation of FSRQs, and FR Is are the parent population of BL 
Lacs (Urry & Padovani 1995). It is believed that the unifi¬ 
cation is mainly due to relativistic beaming. Compared with 
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blazars, radio galaxies is less beaming with a larger jet angle 
to the line of sight. Fanaroff & Riley (1974) separated ra¬ 
dio galaxies into two subclasses with different morphology: 
the peak of low luminosity FR Is is close to the nucleus, 
while high luminosity FR IIs have radio lobes with promi¬ 
nent hot spots and bright outer edges. The distinction of 
morphology also translates into a separation in radio power 
(the fiducial luminosity Lirs Ri 2 X 10^® W Hz-1). ^ 

Owen (1994) found that the FR Is and FR IIs break depends 
on radio and optical luminosity. In addition to host galaxy 
magnitude and environments, FR Is and FR IIs differ in 
optical emission lines (e.g. Zirbel & Baum 1995). As in the 
case of blazars, Ghisellini & Celotti (2001) introduced that 
the division between FR Is and FR IIs actually reflected a 
systematic difference in accretion rate. 

In this paper, after removing beaming effect, we re¬ 
visited the blazar sequence, and made use of the correla¬ 
tion between jet cavity power and iZpeak to study the uni¬ 
fication of blazars and radio galaxies. The paper is struc¬ 
tured as follows: in Sect. 2, we present the samples; the 
results are presented in Sect. 3; discussions and conclu¬ 
sions are presented in Sect. 4. The cosmological parameters 
Hq = 70 km s”! Mpc”!, flm = 0.3 and Ha = 0.7 have been 
adopted in this work. 


2 THE SAMPLES 

Our sample was collected directly from sample of Nemmen et 
al. (2012). But some blazars in Nemmen et al. (2012) were 
not clean Fermi blazars (2LAC; Ackermann et al. 2011). 
We removed the non-clean blazars. Generally, the SED was 
fitted by using a simple third-degree polynomial function. 
However, many blazars were lack of observed SED. Abdo et 
al. ( 2010 b) have conducted a detailed investigation of the 
broadband spectral properties of the y-ray selected blazars 
of the Eermi LAT Bright AGN Sample (LBAS). They assem¬ 
bled high-quality and quasi-simultaneous SED for 48 LBAS 
blazars, and their results had been used to derive empiri¬ 
cal relationships that estimate the position of the two peaks 
from the broadband colors (i.e. the radio to optical, CTro, 
and optical to X-ray, aox, spectral slopes) and from the 7 - 
ray spectral index. Erom 2LACQ, we collected the Qro and 
Qox. Then from the empirical relationships of Abdo et al. 
( 2010 b), we calculated the r/peak and corrected redshift for 
the iZpeak- We excluded the blazars without cTox- At the same 
time, using another empirical formula of Abdo et al. (2010b), 
we estimated the synchrotron peak flux from 5 GHz flux and 
i^peak- The 5 GHz flux is assembled from NASA/IPAG Ex- 
tragalactic Database: NED. When more than one flux was 
found, we took the most recent one. The flux is K-corrected 
according to S,, = ^[(^^(l-l-a)^-^, where a is the spectral in¬ 
dex and a = 0.0. The luminosity is calculated from the rela¬ 
tion = ^TvdiZSv, and du is the luminosity distance. Erom 
Nemmen et al. (2012), for Fermi blazars, /b is the beaming 
factor (/b = 1 — cos(l/r)), where F is the bulk Lorentz fac¬ 
tor of the flow, since blazars obey 9j < l/F (Jorstad et al. 
2005; Pushkarev et al. 2009). Pushkarev et al. (2009) cal¬ 
culated the variability Lorentz factors Pvar from long-term 
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Figure 1. synchrotron peak luminosity versus synchrotron peak 
frequency. Fermi BL Lacs: black empty circles; Fermi FSRQs: 
black filled circles; FR I radio galaxies: red empty triangles; FR 
II radio galaxies: red filled triangles. 

radio observation. The bulk Lorentz factors of Nemmen et 
al. (2012) were collected from the results of Pushkarev et 
al. (2009) (the number ~20%). For blazars without Fvar ob¬ 
tained from Pushkarev et al. (2009), they used the power-law 
fit of /b « 5 X 10 -'*(L4g°)-°'®®^°'^® as an estimator for /(, 
(significant level at the 3.6 ct; is observation 7 -ray lumi¬ 
nosity). Moreover, making use of the relation between cavity 
power and extended radio 300 MHz luminosity and assum¬ 
ing Pjet = Pcav, they estimated the jet kinetic power. So the 
jet power is free of beaming effect. The uncertainty of jet ki¬ 
netic power is 0.7 dex. Our beaming factor (or bulk Lorentz 
factor) and jet kinetic power of blazars are obtained from 
Nemmen et al. (2012). For radio galaxies, we included the 
complete radio galaxies sample of Meyer et al. (2011). The 
r'peak and Lpeak of Meyer et al. (2011) were estimated by fit¬ 
ting non-simultaneous average SED. The jet power of Meyer 
et al. ( 2011 ) was calculated from cavity power as same as 
Nemmen et al. (2012). The relevant data for Fermi blazars 
were listed in Table 1 and for radio galaxies in Table 3 of 
Meyer et al. (2011). 

In total, we have a sample containing 184 clean Fermi 
blazars (98 FSRQs and 86 BL Lacs) and 41 radio galaxies 
(24 FR Is and 17 FR IIs). 


3 THE RESULTS 
3.1 The blazar sequence 

The synchrotron peak luminosity versus synchrotron peak 
frequency is shown in Fig. 1. A distinct “L” or “V” shape 
is not seen in this figure. Blazars and radio galaxies are 
almost located in different regions. The synchrotron peak 
luminosity Apeak of almost all FSRQs (BL Lacs) is larger 
than Apeak of FR IIs (FR Is). From Pearson correlation 
analysis, it is shown that for all blazars sample, there is 
significant anti-correlation while not significant for FSRQs; 
there is low anti-correlation for BL Lacs (see Table 2). In 
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Figure 2. Doppler-corrected synchrotron peak luminosity versus 
Doppler-corrected synchrotron peak frequency for Fermi blazars. 
In the top panel, the data-points are D^-corrected and in the 
bottom panel they are D®-corrected. The meanings of different 
symbols are as same as Fig. 1. 


Fig. 1, there is a BL Lac (J0630.9-2406: logi/peak = 15.39 
Hz, logLpeak = 46.5 erg s“^) which is found as high Lpeak 
and high z^peak blazar in Padovani, Giommi & Rau (2012). 

The Doppler corrections of synchrotron peak frequency 
and luminosity scale as z/ = and Lpeak = Lpeak/S^ re¬ 
spectively, with P = 3 + a for a moving, isotropic source 
and P = 2 + a for a continuous jet (<5 is Doppler factor; 
spectral index a = 1; Urry & Padovani 1995 and Nieppola 
et al. 2008). Since blazars obey 9j < l/P, especially for 
Fermi blazars (Jorstad et al. 2005; Pushkarev et al. 2009; 
Linford et al. 2011; Wu et al. 2014), we assume 5 F (F 
is estimated from radio observation or beaming factor /b). 
From the above relations, Doppler-corrected peak luminos¬ 
ity and frequency can be obtained (D^-correction and D^- 
correction stand for P = 2 -|- a and P = 3 -I- a respectively). 
The Doppler-corrected synchrotron peak luminosity versus 
Doppler-corrected synchrotron peak frequency is shown in 
Fig. 2. The results of Pearson correlation analysis show that 
for all blazar and subclass, there are not significant cor¬ 
relations or are significant positive correlations (both for 
P = 2 + a and P = 3 -I- a; see Table 2). 
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Figure 3. Jet power versus non-beaming-corrected synchrotron 
peak frequency. The uncertainty of jet power corresponds to 
0.7 dex. The discontinuous line corresponds to log Pjet ~ 
44.6 erg s“^. The meanings of different symbols are as same as 
Fig. 1. 


3.2 The unification of blazars and radio galaxies 

Sbarrato, Padovani & Ghisellini (2014) have studied the 
question: “how much do we have to beam the radio lumi¬ 
nosity of the radio galaxies to compare them with blazar?” 
Sbarrato, Padovani & Ghisellini (2014) beamed the radio 
luminosity of the radio galaxies as they were oriented as a 
blazar, who boosted the radio luminosity of the radio galax¬ 
ies by a factor: 

^ 4-/3cosfeL'^ ’ ^ ^ 

where S is the Doppler-boosting factor with apparent su¬ 
perluminal component velocities /3. The Doppler-boosting 
factor is expressed as <5 = with Lorentz factor 

P = (1 — and viewing angle 9. Following the Sbar¬ 

rato, Padovani & Ghisellini (2014), we boost the synchrotron 
peak frequency of radio galaxies by a factor: 

SbL _ 1 — /?COS0RG 
(5rg 1 — /ScosfeL 

We take average viewing angle 6 'rg ~ 40° for radio galax¬ 
ies and 0 bl ~ 3° for blazars similar to Sbarrato, Padovani 
& Ghisellini (2014). If we assume jets of blazar and radio 
galaxy are beamed with a Lorentz factor similar to a com¬ 
mon blazar, then the Lorentz factor P ~ 10. If the jets of 
radio galaxies and blazar are characterized by a rather small 
Lorentz factor, then the Lorentz factor P ~ 3. We also esti¬ 
mate the factor from assuming a average Lorentz factor 
P ~ 5. The factors are 4.7, 11.6 and 37 correspond¬ 
ing to P ~3, 5 and 10 respectively. For radio galaxies, we 
use zZpgg^jj = (|^)i/peak to beam their synchrotron peak fre¬ 
quency (called beam-corrected synchrotron peak frequency). 
The jet power versus non-beam-corrected synchrotron peak 
frequency and beam-corrected synchrotron peak frequency 
can be seen in Fig. 3 and Fig. 4 respectively. From Fig. 3 
and Fig. 4, it is shown that when the synchrotron peak fre- 
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Figure 4. Jet power versus beam-corrected synchrotron peak 
frequency. The beam-corrected synchrotron peak frequency for 
radio galaxy is estimated by assuming P 5 and viewing angle 
^RG ~ 40° for radio galaxies and 0 bl 3° for blazars. The 
error bar is beam-corrected synchrotron peak frequency estimated 
by assuming Lorentz factors of 3 and 10. The uncertainty of jet 
power corresponds to 0.7 dex. The discontinuous line corresponds 
to log ^jet ~ 44.6 erg s The meanings of different symbols are 
as same as Fig. 1. 

quency of radio gala:xy is not beam-corrected, the blazars 
have larger synchrotron peak frequency than radio galaxies. 
Therefore most of blazars do not overlap with radio galax¬ 
ies. However after correcting beaming effect on peak fre¬ 
quency of radio galaxies, we find that most of blazars overlap 
with radio galaxies. From Table 2, we can get that for non- 
beam-corrected peak frequency, there are significant correla¬ 
tions between Tjet and t'peak for blazars-|-radio galaxies, only 
blazars, only BL Lacs, only FSRQs; for beam-corrected peak 
frequency of radio galaxies, there are significant correlations 
between fjet and iZpeak for blazars-l-radio galaxies sample, 
and the correlations of beam-corrected peak frequency are 
much stronger than the correlation of non-beam-corrected 
peak frequency. 


4 DISCUSSIONS AND CONCLUSIONS 

In order to study beaming effect on the blazar sequence and 
obtain jet cavity power, we consider the sample of Nem- 
men et al. ( 2012 ) which is large sample including beam¬ 
ing factor and jet cavity power but not complete sample 
of 2LAC. Finke (2013) studied a 2LAC sample (includ¬ 
ing blazars with unknown redshift), and used the empirical 
equations of Abdo et al. (2010b) to estimate the i/peak and 
Lpeak- Via comparing our Fig. 1 with Fig. 2 of Finke (2013), 
it is found that our Fig. 1 is consistent with Fig. 2 of Finke 
(2013). The correlation analysis of blazars (with known red- 
shift) also supports that there is significant anti-correlation 
between r'peak and Lpeak which is consistent with result of 
Finke (2013). In our Fig. 1, we also include radio galaxies. 
From Meyer et al. (2011), an “L”-shape in the Lpeak — ^^peak 
plot seems to have emerged which destroys the blazar se¬ 


quence. Also Meyer et al. (2011) presented the blazar en¬ 
velope (FR Is and BL Lacs belong to “weak-jet” sources; 
FR IIs and FSRQs belong to “strong-jet” sources). However, 
combining blazars and radio galaxies sample, we hnd that 
our results do not have obvious the blazar envelope, espe¬ 
cially for FR Is+BL Lacs sample. The main reason is that 
there are still blazars in lO^^Hz < iZpeak < lO^^Hz inter¬ 
val for our sample. Our i^peak and Lpeak are estimated from 
the empirical equations of Abdo et al. (2010b) who have 
pointed out that their method assumed that the optical and 
X-ray fluxes are not contaminated by thermal emission from 
the disk or accretion. In blazars where thermal flux compo¬ 
nents are not negligible (this should probably occur more 
frequently in low radio luminosity sources) the method may 
lead to a significant overestimation of the position of Z/'peak- 
From Finke (2013), the following factors can have effect on 
estimating z/peak and Lpeak: the measurement error on de¬ 
termining z.Zpg 3^]5 and Lpeak, the fitting function of SED and 
variability. Finke (2013) have pointed out that the results 
in Lpeak versus z^peak plot are probably not accurate for in¬ 
dividual sources, whereas the overall trend still is present 
even for large errors. From Equations (1), (3) and (4) of 
Abdo et al. (2010b), it is in turn got that a^o oc —logfr, 
logi'peak OC log/r, and log Lpeak OC log/r (/r is radio flux, 
Oro is broadband spectral slopes between radio and optics). 
So even if the correlation between Lpeak and z/peak is dom¬ 
inated due to both dependence on radio flux, the correla¬ 
tion should be positive but not anti-correlation. In addition, 
when accounting for the common dependence on radio flux 
(if any), we use Pearson partial correlation analysis to ana¬ 
lyze the correlation between Lpeak and zzpeak for all blazars. 
The result shows that the significant correlation still exist 
(r = —0.47, P = 3 X lO”*^^). In our sample, there is a BL 
Lac which is found as high Lpeak and high z/peak blazar in 
Padovani, Giommi & Rau (2012). Ghisellini et al. (2012) 
have found that the four high Lpeak and high Vpeak sources 
from Padovani, Giommi & Rau (2012) do not have exces¬ 
sively large 7 -ray luminosity or excessively small 7 -ray pho¬ 
ton index. However, if the high Lpeak and high Z/'peak blazars 
indeed exist, the blazars sequence based on Lpeak — zzpeak 
will be destroyed. 

Another attention of the blazars sequence is beaming 
effect on Lpeak and z.'peak- Then the Doppler factor or bulk 
Lorentz factor should be obtained. The variability Doppler 
factor was derived from the associated variability brightness 
temperature of total radio flux density flares with the intrin¬ 
sic value of brightness temperature, and the bulk Lorentz 
factor from apparent super-luminal component velocities 
and variability Doppler factor (Hovatta et al. 2009). For the 
blazars without without Fvar obtained from Pushkarev et al. 
(2009), Nemmen et al. (2012) used 7 -ray luminosity to esti¬ 
mate beaming factor with 3.6 (t significant level. For blazars 
with variability bulk Lorentz factor available, the relative 
uncertainty in F is 0.3 and the relative uncertainty in P es¬ 
timated from the Lj — /b relationship is 0.7 (Nemmen et al. 
2012). Furthermore, we assume F ~ <5 which is reasonable 
because 6j < l/F for blazars, especially for Fermi blazars 
with a larger beaming factor and smaller Oj (Jorstad et al. 
2005; Pushkarev et al. 2009; Linford et al. 2011; Wu et al. 
2014). It is noted that the above estimating S is probably 
not accurate for individual sources. From Fig. 2 and Table 
2, we hnd that after being Doppler-corrected, for all blazar 
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and subclass, the correlations between Lpeak and i^peak are 
not significant or become significant positive correlations, 
i.e. the anti-correction between Lpeak and t'peak disappears, 
which is consistent with result of Nieppola et al. (2008). 

Through studying the correlation between jet power 
and the synchrotron peak frequency, it is found that there is 
significant anti-correlation between -Fjet and Z/'peak for Fermi 
blazars, which supports the blazar sequence, i.e. a stronger 
radiative cooling for higher jet power sources results in 
smaller energies of the electrons emitting at the peaks. In 
addition, the results also imply the BL Lacs - FSRQs di¬ 
vide. From Ghisellini et al. (2009, 2011, 2015), Sbarrato et al. 
(2012), Sbarrato, Padovani & Ghisellini (2014) and Xiong & 
Zhang (2014), it is got that Pjet ~ Ldisk and the physical dis¬ 
tinction between BL Lacs and FSRQs is LeLu/LEdd ~ 10“^. 
Based on result of Xiong & Zhang (2014), the average black 
hole mass of Fermi blazars is ~ lO^ ^M©. Then the jet power 
of BL Lacs and FSRQs divide corresponds to log Pjet ~ 
44.6 erg s“^(LEdd = 1.3x 10 ®®( j||-)erg s“\ Ldisk « IOLblr, 
M ^ 10® ®Mq). The Fig. 3 and Fig. 4 present the divide line 
which also supports the physical distinction. The BL Lacs 
with logPjet > 44.6 erg s“^ (i.e. LeLn/LEdd > 10“®) can be 
considered as transition sources (see Ghisellini et al. 2011). 
Moreover, the BL Lac which is found as high Lpeak and 
high iZpeak blazar in Padovani, Giommi & Rau (2012) does 
not have high jet cavity power (logPjet = 43.2 erg s“^). A 
possible explanation is that the high Lpeak and high zzpeak 
blazar is resulted from beaming effect while the blazar virtu¬ 
ally is low power source. Ghisellini et al. (2009) have studied 
the 7 -ray luminosity versus y-ray photon index (a proxy for 
inverse Gompton peak frequency), and proposed the Fermi 
blazars divide. We should note that there is beaming effect 
on 7 -ray luminosity while not on jet power estimated from 
extended radio emission. 

Next, we will discuss the unification of blazars and ra¬ 
dio galaxies. The Fig. 3 and Fig. 4 show that when the 
synchrotron peak frequency of radio galaxy is not beam- 
corrected, most of blazars do not overlap with radio galax¬ 
ies (especially for FR Is and BL Lacs), while after being 
beam-corrected peak frequency for radio galaxies, most of 
blazars overlap with radio galaxies. Nieppola et al. (2008) 
have found the dependence between Doppler factor and 
r'peak for blazars. So the Vpeak should be Doppler-corrected 
or ruled out the influence of beaming effect. However, it is 
not possible to derive the specific Lorentz factor/Doppler 
factor for each source, especially for radio galaxies. In this 
paper, we adopt the method of Sbarrato, Padovani & Ghis¬ 
ellini (2014). Making use of different Lorentz factor, we ana¬ 
lyze the correlations between jet power and beam-corrected 
peak frequency for blazars-fradio galaxies sample. From Ta¬ 
ble 2 , we can get that after being beam-corrected peak fre¬ 
quency of radio galaxies, there are significant correlations 
between jet power and peak frequency for blazars-fradio 
galaxies sample, and the correlations of beam-corrected peak 
frequency are much stronger than the correlation of non- 
beam-corrected peak frequency. When beaming-corrected 
peak frequency, the assuming Lorentz factors are reason¬ 
able. For blazars, the F ~3, 5 and 10 are possible (see Sbar¬ 
rato, Padovani & Ghisellini (2014) for detail). Giovannini 
et al. ( 2001 ) studied a complete sample of radio galaxies, 
and estimated that relativistic jets with Lorentz factor in 
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Figure 5. 7 -ray photon index versus jet power. The uncertainty 
of jet power corresponds to 0.7 dex. The meanings of different 
symbols are as same as Fig. 1. 

the range 3-10 are present in high and low power ra¬ 
dio sources. Now, we explain how to associate the correla¬ 
tions between jet power and beam-corrected r'peak with the 
unification of blazars and radio galaxies. Urry & Padovani 
(1995) have explained that the radio extended luminosity 
of FR IIs (FR Is) and FSRQs (BL Lacs) should be com¬ 
parable. Our jet cavity power is estimated from radio ex¬ 
tended luminosity. Then the jet power between FR IIs (FR 
Is) and FSRQs (BL Lacs) should be comparable. In addi¬ 
tion, as discussed in the previous paragraph, the division 
between FSRQs and BL Lacs is log Pjet ~ 44.6 erg s ^ or 
Ld/LEdd ~ 10“^. According to the two-dimensional optical- 
radio luminosity correlation from Owen & Ledlow (1994), 
Ghisellini & Celotti (2001) found that the division between 
FR IIs and FR Is is Pjet ~ 0.015LEdd. If the radio galax¬ 
ies have similar average black hole mass with blazars, then 
the division between FR IIs and FR Is also is log Pjet ~ 
44.6 erg s“^. From Fig. 3 and Fig. 4, it is seen that the line 
of logPjet ~ 44.6 erg s“^ also can separate FR Is from FR 
IIs. Abdo et al. (2010b) subdivided blazars based on the 
I'peak (low/intermediate/high-synchrotron-peaked blazar). 
Almost all FSRQs are low-synchrotron-peaked blazar and 
BL Lacs include low/intermediate/high-synchrotron-peaked 
blazar. Based on unified scheme, after being beam-corrected, 
the iZpeak should be comparable between FR IIs (FR Is) and 
FSRQs (BL Lacs) because the differences between FR IIs 
(FR Is) and FSRQs (BL Lacs) mainly rely on jet angle to 
the line of sight. Therefore based on above discussion, in 
Pjet — r'peak plot, blazars and radio galaxies should have a 
similar correlation, i.e. the blazars and radio galaxies can be 
unified from the correlation between jet power and beam- 
corrected peak frequency. At present, since only there are 
limited radio galaxies with precise Lorentz factor estimated 
from observation, it is hard to accurately estimate beam- 
corrected peak frequency for individual sources. In order to 
test whether assuming Lorentz factors bring in remarkable 
effect on the overall trend which results in spurious corre¬ 
lation, We also analyze the correlation between jet power 
and 7 -ray photon index for Fermi blazar and radio galaxie 
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because 7 -ray photon index is free of beaming effect and 
Lj versus 7 -ray photon index has been used to support the 
blazar sequence. If for blazars and radio galaxies, the corre¬ 
lations between jet power and 7 -ray photon index still exist 
and are consist with the correlation between Pjet and beam- 
corrected i^peaki then it is clarihed that our assuming Lorentz 
factors do not bring in remarkable effect on the overall trend. 
From Fig. 5 (FR Is: 3C 78, 3C 84, M87, Cen A, NGC 6251; 
FR IIs: 3C 380, 3C 111), the correlations between jet power 
and 7 -ray photon index still exist and support the blazar 
sequence (higher jet power - softer photon index). In the 7 - 
ray photon index versus jet power plot, FR IIs overlap with 
FSRQs and most of FR Is with BL Lacs, which support the 
unification of blazars and radio galaxies, similar to the cor¬ 
relation between Pjet and beam-corrected r'peak- Therefore, 
the assuming Lorentz factors do not bring in remarkable 
effect on the overall trend. 
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Table 1. The Fermi blazar sample. 


2FGL name 

Type 




-^R 

1 (®) 


J0050.6-0929 

BLL 

0.634 

44.31 

-2.08 

1920 

14.83 

46.86 

J0108.6+0135 

FSRQ 

2.099 

46.46 

-3.21* 

3400 

13.10 

47.08 

J0112.1+2245 

BLL 

0.265 

43.75 

-2.27 

327 

14.98 

45.52 

J0112.8+3208 

FSRQ 

0.603 

45.21 

-2.52 

405 

13.15 

45.37 

J0116.0-1134 

FSRQ 

0.67 

45.46 

-2.43 

1900 

13.30 

46.13 

J0120.4-2700 

BLL 

0.559 

45.06 

-2.45 

1035 

13.99 

46.10 

J0132.8-1654 

FSRQ 

1.02 

45.22 

-2.65 

1584 

13.73 

46.60 

J0136.9+4751 

FSRQ 

0.859 

44.78 

-2.62* 

3150 

13.69 

46.73 

J0137.6-2430 

FSRQ 

0.838 

45.57 

-2.45 

1558 

13.58 

46.37 

J0141.5-0928 

BLL 

0.733 

45.08 

-2.41 

1200 

14.06 

46.41 

J0205.3-1657 

FSRQ 

1.739 

45.42 

-2.79 

1400 

13.46 

46.79 


(a) redshift directly collected from NED. (b) jet power estimated by cavity power in unit of erg s“^. (c) beaming factor /b = 1 —cos(l/r). 
values with a “ * ” represent that they are directly from radio observation, (d) total radio 5 GHz flux in mjy. (e) rest-frame synchrotron 
peak frequency and luminosity in Hz and erg s“^. 

(f) This table is published in its entirety in the electronic edition. A portion is shown here for guidance. 


Table 2. The results of correlation analysis. 


Sample 

-^peak 

(a) 

versus \ 

peak 

r D2 D(b) 

L , versus i/ ^ / 

peak peak 

T-D3 t)(c) 

L ^ , versus u ^ / 

peak peak 

D (®) 

Pjet versus 


r(f) 

p(f) 

r P 

r P 

r 

P 

Blazars 

-0.5 

4.2 X 10-13 

0.009 0.9 

0.26 3 X 10-1 

-0.69 

3 X 10-31 

BL Lacs 

-0.24 

0.02 

-0.03 0.8 

0.12 0.26 

-0.45 

1.7 X 10-3 

FSRQs 

-0.18 

0.07 

0.33 1 X 10-3 

0.46 1.5 X 10-s 

-0.52 

4.8 X 10-3 

Sample 

Pjet versus Ffffl (F ~ 3) 

Pjet versus ~ 5) 

Pjet versus (F ~ 10 ) 

D (®) 

Pjet versus 

Blazars-hRG 

-0.59 

2.6 X 10“33 

-0.62 1.2 X 10-33 

-0.64 2.1 X 10-31 

-0.50 

2.6 X 10-13 


(a) Synchrotron peak luminosity versus synchrotron peak frequency, (b) Doppler-corrected synchrotron peak luminosity versus Doppler- 
corrected synchrotron peak frequency (Z)^-corrected). (c) Doppler-corrected synchrotron peak luminosity versus Doppler-corrected syn¬ 
chrotron peak frequency (D^-corrected). (d) jet power versus beam-corrected synchrotron peak frequency from assuming F ~3, 5 and 
10 . (e) jet power versus non-beam-corrected synchrotron peak frequency. 

(f) The r is the Pearson correlation coefficient; P is the chance probability. 








